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Abstract. Magnetic and charge susceptibilities of the two-dimensional repulsive Hubbard 
model are investigated applying a strong coupling diagram technique in which the expansion 
in powers of the hopping constants is used. For small lattices and high temperatures results 
are in agreement with Monte Carlo simulations. With the departure from half-filling x the 
low-frequency magnetic susceptibility becomes incommensurate and the incommensurability 
parameter grows with x. The incommensurability, its dependence on frequency and on x 
resemble experimental results in lanthanum cuprates. Also for finite x sharp maxima appear 
in the static charge susceptibility. The maxima are finite which points to the absence of the 
long-range charge ordering (static stripes). However, for x w 0.12 the maxima are located 
near the momenta (0, ±tt/2), (±7r/2, 0). In this case an interaction of carriers with tetragonal 
distortions can stabilize stripes with the wavelength of four lattice spacings, as observed in 
the low-temperature tetragonal phase of cuprates. As follows from the obtained results, the 
magnetic incommensurability is not a consequence of the stripes. 



The Hubbard model is thought to be appropriate to describe the main features of electron 
correlations in narrow energy bands, leading to collective effects such as magnetism and 
metal-insulator transition. It has often been used to describe real materials exhibiting these 
phenomena. In more than one dimension, the model is not exactly solvable and a variety of 
numerical and analytical approximate methods was used for its study. In the case of strong 
electron correlations inherent in cuprate perovskites, when the ratio of the hopping constants 
to the Hubbard repulsion is a small parameter, the strong-coupling diagram technique [TJ [2] 
can be applied for the investigation of the model. In this technique, Green's functions are 
calculated using expansions in powers of the hopping constants. The terms of these expansions 
are expressed by means of site cumulants of electron creation and annihilation operators. 

In this article we use this approach for calculating magnetic and charge susceptibilities of the 
2D repulsive Hubbard model which presumably describes the Cu-0 plane of cuprate perovskites. 
Our aim is to investigate possible mechanisms of the incommensurate magnetic response [3] and 
the stripe formation [3j observed in these crystals. The susceptibilities are connected with the 
spin and charge Green's functions 

D{\'t,\t) = {Ts^{t')s^{t)), B(l'r, It) = (T5n v {T')5n x {T)), (1) 



where the angular brackets denote the statistical averaging with the Hubbard hamiltonian which 
also determines the time evolution of operators in Eq. (pQ), T is the chronological operator, the 



spin operators = a] ai- a , a\ a and a\ a are the electron creation and annihilation operators on 
the site 1 of the square lattice with the spin projection a, 8n\ = n\ — (n\) with n\ = a \ a a ia- 
Using the strong-coupling diagram technique one can convince oneself that equations for D 
and B look similar jjj]. For the function B they read 
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Here the combined indices q = (k, iui v ) and pj = (]Cj,iu> n .) were introduced, u v = 2wkT and 
oo n = (2n + 1)ttT are the boson and fermion Matsubara frequencies with the temperature T, k 
is the wave vector, G(p) = ((akaW^}) is the electron Green's function, Tl(p) = 1 + tkG(p), is 
the Fourier transform of the hopping constants, Q(p) = t\JJ(p) is the renormalized hopping, 
N is the number of sites, A(pi,p + pi,p + P21P2) is the sum of all four-leg diagrams and 
^{PiiP + Pi,P + P21P2) is its irreducible subset. With explicitly specified spin indices of the 
external lines, this latter notation reads A = \ J2a'cr M "' \ G ' 1 a i a )- The difference between Eq. ([2]) 
and the equations for D is in the spin indices of this quantity. In the latter case A(|, j, j, |) 
enters into the equations. 

In the calculations we approximated quantities II(p) and <d(p) by 1 and tu_, respectively, 
and used the lowest-order irreducible four-leg diagrams - the second-order cumulants - for A 
in Eq. ([2|). A comparison of the obtained results with Monte Carlo simulations is shown in 
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Figure 1. The static charge susceptibility in a 
4x4 lattice for t = -U/8 and T = 0.125|i|. The 
susceptibility is plotted along the triangular 
contour in the Brillouin zone. The corners of 
the contour are given by the momenta k = 
(0,0) (r), (vr,0) (X), and (vr,7r) (M). Open 
symbols are results of Monte Carlo simulations 
for n = 1 (squares) and n = 0.95 (circles). 
Filled symbols are our results for the same 
electron fillings. 



Fig. [TJ In this figure, t is the nearest-neighbor hopping constant, U is the Hubbard repulsion 
and n = {n{) is the electron filling. 

For small real frequencies uj the momentum dependence of the imaginary part of the magnetic 
susceptibility x"(kuj) = ImD(kj) calculated in a 100x100 lattice is shown in Fig. [2j At half- 
filling, n = 1, for the considered temperature T = 0.05U the system is in the state with 
the commensurate short-range order [the dash-dotted curve in Fig. [2](b)] . With departure 
from half-filling, the maximum of ~x"(\z.oj) shifts from (tt, tt) and the susceptibility becomes 
incommensurate. For x = |1 — n\ > 0.06 the susceptibility is peaked at momenta (vr, 7r ± 5), 
(it ± 5, it) [see Fig. E|a)] and the incommensurability parameter d grows with x [see Fig. [2](b)] . 
For a fixed n, 5 decreases with increasing uj and, at some frequency u r , the incommensurability 
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Figure 2. (a) The momentum dependence of x"(kw) along the edge [solid line, k = (rr, n)] 
and diagonal [dashed line, k = (n, k)] of the Brillouin zone for t = — O.llf/, u = 0.002^7 and 
n ~ 0.88). (b) The momentum dependence of x"(kw) along the zone edge for n ~ 0.88 (solid 
line), n rs 0.94 (dashed line), and n = 1 (dash-dotted line), t = -0.11U and to = 0.002f/. (c) 
The dispersion of maxima in x"(ku;) along the zone edge for t = —0.11U and n ~ 0.88. 



disappears and x"(kw) appears to be peaked at (tt, it) [see Fig. EJ^c)]. The same behavior of the 
magnetic susceptibility at low frequencies was observed in lanthanum cuprates [3|- In Fig. [21 the 
values of 5 are close to those measured experimentally for the same x. Parameters of this figure 
were chosen so that for n ~ 0.88 the value of u r is close to the experimental value 50 meV. For 
the superexchange constant J = 4t 2 /U = 0.15 eV we find ui r = 44 meV from Fig. [2j 

The momentum dependence of the static charge susceptibility Xc(k) = -B(k, v = 0) is shown 
in Fig. [3] for two electron fillings. At half-filling the susceptibility is small and its dependence 
on momentum is weak. Immediately after the Fermi level crosses one of the Hubbard subbands, 
which leads to departure from half-filling, a sharp peak appears in x c (k) near the T point [see 
Fig. E^a)]. With increasing x the susceptibility grows and the peak transforms to a ridge around 
r [see Fig. [3](b)] . Notice that in the entire considered range < x < 0.2 the static charge 
susceptibility remains finite. This means that a long-range or stripe ordering of charges does not 
occur in the Hubbard model for the considered range of parameters corresponding to cuprate 
perovskites. The variation of the ratio U/t in this range and the inclusion of the hopping to 
more distant sites do not change this conclusion. 

This result indicates that the above-discussed magnetic incommensurability is not connected 
with charge stripes and is a consequence of strong electron correlations. 

Although the above results do not support a purely electronic mechanism of long-range charge 
ordering, they give some insight into the way in which phonons can stabilize stripes. An essential 
role of certain CuOg octahedra tilts in such stabilization follows from the fact that static stripes 
were observed only in the low-temperature tetragonal (LTT) phase of lanthanum cuprates [1]. 
It is known that such tilts are strongly coupled to the carriers [7]. This interaction leads to the 
following term in the adiabatic potential: 

A£ = -V Xc (k)^, (3) 
z k 



where the function 0k contains only even powers of tilt coordinates. As follows from Eq. (|3j), 



Figure 3. The static charge susceptibility in a 40x40 lattice for t = —U/8 and T = 0.001C/. 
The susceptibility is shown in the first quadrant of the Brillouin zone for the electron fillings 
n = 0.97 (a) and 0.88 (b). 



finite values of the vibration coordinates give an energy gain and the larger the value of the 
susceptibility is, the larger energy gain can be achieved. As seen from Fig. [HJ^b), for x ~ 0.12 
the susceptibility maxima are located on the axes of the Brillouin zone approximately halfway 
between its center and its boundary. Such maxima give the lowest energy gain for the charge 
density wave with the wavelength equal to four lattice spacings, as observed in the LTT phase 
of lanthanum cuprates |3]. In contrast to other electron fillings such a wave is commensurate 
with the lattice which is essential for its stability. 

In summary, the magnetic and charge susceptibilities of the 2D repulsive Hubbard model 
were calculated using the strong-coupling diagram technique. It was found that with 
departure from half-filling the low-frequency magnetic susceptibility becomes incommensurate. 
This incommensurability, its dependence on frequency and on electron filling are similar to 
experimental results in lanthanum cuprates. In the range of parameters corresponding to cuprate 
perovskites the maxima of the static charge susceptibility are finite which points to the absence 
of the long-range charge ordering (static stripes) in the Hubbard model. However, the shape of 
the obtained susceptibility for x ~ 0.12 suggests that an interaction of carriers with tetragonal 
distortions can stabilize stripes with the wavelength of four lattice spacings, as observed in the 
low-temperature tetragonal phase of cuprates. From the obtained results it follows that the 
magnetic incommensurability is not a consequence of the stripes. 
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